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We present the results of the first exposure of a Liquid Argon TPC to a multi-GeV neutrino beam. 
The data have been collected with a 50 liters ICARUS-like chamber located between the CHORUS 
and NOMAD experiments at the CERN West Area Neutrino Facility (WANF). We discuss both 
the instrumental performance of the detector and its capability to identify and reconstruct low 
multiplicity neutrino interactions. 

PACS numbers: 29.40.Gx, 13.15.+g 



I. INTRODUCTION 



The Liquid Argon Time Projection Chamber (LAr 
TPC) technology was originally proposed by C. Rubbia 
in 1977 0, as a novel concept for a massive neutrino 
detector that allows three-dimensional event reconstruc- 
tion with high accuracy and precise calorimetric mea- 
surement. The feasibility of the technology, developed 
inside the context of the ICARUS project Q, has been 
demonstrated by an extensive R&D programme that ad- 
dressed the main technological challenges of the LAr de- 
tection technique (i.e. argon purification, wire chambers, 
electronics, DAQ, etc.). The ICARUS collaboration has 
been able to operate prototypes of increasing mass (3 
tons 0, 0|, 14 tons 0], 600 tons 0) and therefore the 
LAr TPC can now be considered as a mature technology 
able to produce multi-kton detectors for astroparticle and 
neutrino physics. 



Within the broad physics programme of the ICARUS 
project, an important issue is the study of oscillations 
through the detection of long-baseline neutrino interac- 



tions. The possibility to observe v„ 



and v., 



by means of kinematic criteria is known to be limited, 
among others, by the knowledge we have of nuclear ef- 
fects (Fermi motion, nuclear re-scattering and absorp- 
tion, etc.). Therefore it is very important to acquire 
experimental data in order to tune the existing Monte- 
Carlo models. In addition, the exposure of a LAr TPC to 
a neutrino beam is mandatory to demonstrate the high 
recognition capability of the technique and gain experi- 
ence with real neutrino events. 

In 1997, the ICARUS collaboration together with a 
group from INFN and Milano University Q proposed to 
expose a 50 liter LAr TPC to the multi-GeV wide band 
neutrino beam of the CERN West Area Neutrino Facility 
(WANF) 0, during the NOMAD and CHORUS pjj 
data taking. The test was part of an R&D program for 
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a medium baseline v r appearance experiment 11]. The 
idea was to collect a substantial sample of quasi-elastic 
interactions [y^ + n — > p + fi~) to study the following 
physics items: 

• Measurement of the acoplanarity and missing 
transverse momentum in events with the fj,-p topol- 
ogy in the final state, in order to assess Fermi mo- 
tion and proton re-scattering inside the nucleus. 

• Appearance of nuclear fragments (short tracks and 
blobs around the primary interaction vertex) in 
quasi-elastic events. 

• A preliminary evaluation of e/7 and e/ir° discrimi- 
nation capability by means of the specific ionization 
measured on the wires at the beginning of the can- 
didate track. This measurement is limited by the 
size of the chamber. 

The data collected in 1997 offer the unique opportu- 
nity to study nuclear effects in the Argon nucleus and 
to assess the identification and reconstruction capability 
of a LAr TPC for low- multiplicity neutrino events. The 
present work shows, for the first time, a comprehensive 
set of results from the 1997 test (see [!l E 111 E3 for 
preliminary results). The structure of the paper is as fol- 
lows: the experimental setup and the instrumental per- 
formance of the 50 liters TPC are discussed in Sec. HH and 
Sec. IIIII respectively. In Sec. lIVI the event reconstruction 
and particle identification are detailed. The analysis of a 
"golden-sample" of quasi-elastic (QE) CC interactions 
is presented in Sec. together with a comparison with 
theoretical expectations. 

II. THE EXPERIMENTAL SETUP 

The LAr TPC was placed on a platform 4.5 meters 
above ground, right in between the CHORUS and NO- 
MAD detectors (Fig. [TJ). Tne modest size of the LAr 
TPC fiducial volume (^50 liters), coupled with the high 
energy of the WANF v beam, made necessary a muon 
spectrometer downstream the TPC. A coincidence with 
the NOMAD DAQ was set up to use the detectors located 
into the NOMAD magnetic dipole as a magnetic spec- 
trometer. The experimental setup was completed with 
additional counters for the trigger and veto systems. 

The LAr TPC (Fig0) has an active volume of 
32x32x46.8 cm 3 , enclosed in a stainless steel vessel in 
the shape of a bowed-bottom cylinder 90 cm high and 
with a radius of 35 cm. The active volume corresponds 
to 67 kg of Liquid Argon (T = 87K at 1 atm, p = 1.395 
g/cm 3 ). Ionization electrons produced by the passage 
of charged particles drift vertically toward the anode by 
means of a constant electric field of 214 V/cm. The read- 
out area (anode) is made up of two orthogonal wire planes 
mounted at a distance of 4 mm. Each stainless steel 
wire has a diameter of 100 //m; the distance between the 
wires is 2.54 mm and each plane is made up of 128 wires. 



The cathode and the lateral field-shaping electrodes are 
copper strips (5 mm thickness) positioned on a vetron- 
ite support with printed board techniques. The support 
was glued on a honeycomb structure to ensure rigid- 
ity. The distance between two adjacent strips is 10 mm. 
The 10 KV drift voltage is distributed to the strips by 
means of 100 Mil resistors. The ionization electrons drift 
through the first wire plane (induction plane) so that the 
integrated signal induced in it is zero. The second plane 
(collection plane) stops the electrons and acts as charge 
collector. The mean charge per unit pitch (2.54 mm) for 
a mip crossing the chamber horizontally corresponds to 
about 2.3 x 10 4 electrons (~4 fC) assuming no charge loss 
along the drift volume. The smallness of the signal and 
the absence of an amplification phase requires very low 
noise charge preamplifiers. They were based on JFET 
transistors and have been operated into the liquefied gas 
in the proximity of the readout wires. The polarization 
voltage (15 V) is distributed separately to each preampli- 
fier to minimize the number of dead channels in case of 
failure of one of the components during data taking. The 
output signal is transmitted up to the amplifiers, which 
have been operated in current mode and are located at 
room temperature. The amplified signal is brought up to 
a set of fast 8-bit ADC. This configuration [lj] allowed 
a signal to noise ratio of 11 with mip signals equivalent 
to 10 ADC counts and no saturation even in the occur- 
rence of e.m. showers. The signal is sampled with a 
2.5 MHz frequency for a duration of 500 fis, correspond- 
ing to the highest possible drift time (primary ionization 
at the cathode). The stream is recorded into a buffer. 
The arrival of a subsequent trigger causes the switch of 
the data stream to another buffer (up to 8 buffers are 
available). The multi-buffer writing procedure minimize 
the dead time for signal recording to less than one sam- 
pling time (i.e. 500 (is). 

As mentioned above, the active part of the detector 
is located inside a stainless steel cylinder. The connec- 
tion to the outside area is obtained through a set of UHV 
flanges housing the signal, the high voltage cables and the 
vacuum feed-through. The cylinder is positioned into a 
1 m diameter dewar partially filled with low purity Liq- 
uid Argon acting as a thermal bath and it is rotated 30° 
along the vertical axis with respect to the nominal beam 
direction to reduce the number of particles crossing just 
one readout wire. The ceiling of the external dewar is 
in direct contact with air. The level of Argon in the 
dewar is so small that the pure Ar in the inner cylin- 
der can evaporate as well. The gaseous Ar in the active 
volume crosses the feed-through and reaches an Oxisorb 
filter for purification [4j. After, it is liquefied inside a 
second buffer cooled by low-purity Ar and, finally, recol- 
lected into the active volume. During this test, the Argon 
has been doped with Tetramethyl Germanium in order 
to partially reconvert scintillation light into ionization 
flijj ]. Continuous recirculation of the Argon through the 
purifier keeps the level of contamination stable against 
micro-leaks or outgassing. The total Argon consumption 
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FIG. 1: A sketch of the experimental setup (top view). Relative sizes are not to scale. 




FIG. 2: The 50 liters Liquid Argon Time Projection Chamber. 



necessary to circulate the pure liquefied gas and to com- 
pensate for the heat losses was ^200 liters per day. The 
performance of the recirculation system is described in 
Sec. EH 

The chamber has been exposed to the v beam produced 
at the CERN West Area Neutrino Facility (WANF H 
E3). The primary 450 GeV protons from the CERN SPS 
are extracted every 14.4 s in two spills of 6 ms duration 
each and separated by 2.5 s. On average, 1.8 x 10 13 pro- 
tons per spill hit a segmented beryllium target. Secon- 



daries are selected in momentum and focused by a system 
of collimators and magnetic lenses. They reach a 289.9 m 
decay tunnel followed by an iron and earth absorber. The 
experimental area is located 835 m downstream the tar- 
get. The mean energy of the v ^ reaching the detectors is 
24.3 GeV when integrating over the active NOMAD area 
(2.6 x 2.6 m 2 ), while contaminations from other flavors 
are below 7% for and ~ 1% for v e 0]. Being the 
TPC centered on the beam axis and covering a smaller 
surface, an harder neutrino flux is expected (with a mean 
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energy of about 30 GeV). 

The trigger is provided by a set of scintillators lo- 
cated between the chamber and the NOMAD apparatus 
(Fig. . Each of the 3 trigger scintillator counters has a 
dimension of 110x27 cm 2 . They are positioned in a half- 
circle 60 cm beyond the center of the chamber. Incoming 
charged particle are vetoed by 5 scintillators mounted in 
front of the chamber, 50 cm before its center, and by the 
last scintillator plane of CHORUS. The latter vetoes par- 
ticles deflected by the CHORUS magnetic field entering 
the chamber at large angles with respect to the nominal 
beam direction. 

The trigger requires the coincidence of the SPS beam 
spill and at least one of the trigger scintillators, vetoed 
by the scintillators put in front of the chamber and the 
CHORUS plane. This trigger is put in coincidence with 
the two trigger scintillator planes of NOMAD (Tl and T2 
in Fig.QJ [l||. Moreover, a trigger is rejected if the NO- 
MAD acquisition system is in BUSY mode or if the delay 
with respect to the previous trigger is lower than 500 [is 
("drift protection") 1 . The kinematic reconstruction of 
the outgoing muon from NOMAD (see Sec llV Djl allowed 
us to identify the events produced outside the fiducial 
volume that trigger the chamber. The vast majority are 
v interactions in the dewar or in the thermal bath plus 
a contamination of crossing muons due to veto inefficien- 
cies. Note that the request of a charged particle trigger- 
ing the chamber locally and reaching NOMAD up to Tl 
and T2 inhibits the acquisitions of neutral-current and v e 
charged-current events, limiting the sample of the present 
test to v^i charged-current interactions. The trigger effi- 
ciency has been monitored during data taking through a 
dedicated sample of through-going muons. It turned out 
to be 97% averaged over the whole data taking period. 



at the level of 0.1 ppb to allow drifts over O(oo) meters. 
The electron lifetime r e is defined by: 



Q(t) = Qoe"^ (1) 

Qo being the primary deposited ionization charge and 
t the electron drift time. The to is provided by the trigger 
scintillation counters. 

The effectiveness of the purification system through 
direct measurements of r e has been monitored during 
data taking by a dedicated setup first developed by the 
ICARUS Collaboration in 1989 0. It consists of a small 
double-gridded drift chamber located below the readout 
planes. Electrons are generated near the cathode via pho- 
toelectric effect driven by a 20 ns UV laser pulse. Each 
pulse produces a bunch of 10 7 electrons. They drift to- 
ward the anode along the electric field lines and cross a 50 
cm drift region between the two transparent grids. The 
ratio of the induced signal in the proximity of the grids 
provides a real time estimate of r e during data taking 
and, hence, a measurement of the Ar contamination (see 
Fig. |3J . An independent estimate of r e can be obtained 
monitoring the charge attenuation for muons crossing the 
chamber at various heights. Both methods provided con- 
sistent results and the electron lifetime averaged along 
the whole data taking period turned out to be always 
higher than 10 ms. This performance was excellent since 
the maximum drift time was about 400 [is and therefore 
the attenuation of the ionization over the drift distance 
was negligible. 



III. INSTRUMENTAL PERFORMANCE OF 
THE 50 LITERS TPC 

A. Liquid Argon purification and electron lifetime 

A crucial working parameter for a TPC operating with 
liquefied rare gases is the lifetime of free electrons in the 
medium. Primary electrons produced by the passage of 
charged particles drift toward the anode crossing macro- 
scopic distances. In the present case, for an ionization 
electron created in the proximity of the cathode, the drift 
path-length exceeds 46 cm. The effectiveness of charge 
collection at the readout plane is related to the purity 
of the Argon since electron-ion recombination is m ainly 
due to oxygen molecules present in the LAr bulk |19| . 
The contamination of electronegative molecules must be 



This condition inhibits the occurrence of event overlaps in the 
multibufFer readout system. 



B. Drift velocity 

An absolute determination of the drift velocity and 
its uniformity along the fiducial volume of the chamber 
has been obtained exploiting the external scintillators 
and the additional information from NOMAD. A dedi- 
cated trigger selecting through-going muons ( "imp sam- 
ple") has been put into operation adding two additional 
scintillators (30x30 cm 2 ) before and after the external 
dewar. The relative position of the chamber with respect 
to the NOMAD reference frame has been obtained by 
residual minimization of the track parameters recorded 
both by the TPC and the NOMAD data acquisition sys- 
tem [lj|- The drift velocity is obtained by fitting the 
absolute vertical position of the muon as reconstructed 
by NOMAD trackers versus the drift time measured lo- 
cally by the TPC. The corresponding drift velocity is 
Vd = 0.905 ± 0.005 mm//xs (this value differs from the 
one quoted in [2lJ due to different running conditions of 
pressure and temperature). No systematic biases have 
been observed in the fiducial volume used for the neu- 
trino interaction analysis (Sec[Vj). 
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FIG. 3: Lifetime of the drifting electrons in the 50 liters Liquid Argon TPC at the CERN v beam. The filling of the chamber 
with LAr was performed on the fourth of April 1997. 



IV. EVENT RECONSTRUCTION AND 
PARTICLE IDENTIFICATION 



A. Spatial reconstruction 

Interactions happening in the TPC fiducial volume 
were fully registered in two 2D images with a common 
coordinate (time), with full calorimetric information as- 
sociated to each point. Each 2D-image represents the 
signal amplitude digitized by the ADC (in a linear gray 
scale) versus the time sample (drift coordinate) and the 
wire number. Fig. 0] and Fig. [5] (top part) show the raw 
images of two CC events in the collection (left panel) 
and induction plane (right panel). These raw images rep- 
resent the time evolution (in vertical axis) of the signal 
induced in the wires (in horizontal axis). 

The purpose of the reconstruction procedure is to ex- 
tract the physical information provided by the wire out- 
put signals, i.e. the energy deposited by the different 
particles and the point where such a deposition has oc- 
curred, to build a complete 3D and calorimetric picture 
of the event. In this sense, the spatial reconstruction 
of different tracks is needed first in order to compute 
the calibration factors entering in the calorimetric recon- 
struction of the events, as explained in Sec. IIV Bl A de- 
tailed description of the spatial reconstruction tools was 
reported elsewhere 0, ■ For the present analyses new 
significant improvements have been carried out |23| . 

The basic building block of a track is called a hit, de- 
fined as the segment of track whose energy is detected by 
a given wire. A precise determination of the hit position 
and charge is carried out fitting the wire signal to an ana- 
lytical function which describes the detector response . 
Identified nearby hits are associated into 2-dimensional 
clusters. Track candidates are searched for by means of a 
tree algorithm as in p3 | and/or using a neural network- 



based algorithm as in j25j. Both methods efficiently de- 
tect interaction vertexes and finally, split clusters into 
smooth 2D tracks. Finally, 3D tracks are built matching 
2D track projections from both views. 

In Fig. 21 we show an example of the full reconstruc- 
tion procedure for a low- multiplicity CC event. This 
algorithm allows the reconstruction of more complicated 
event topologies (see Fig. EJ). This is the first time that 
interactions of multi-GeV accelerator neutrinos occurring 
in a LAr TPC are fully reconstructed in 3D. 



B. Calorimetric reconstruction 

The ionization charge is precisely measured at the col- 
lection wire plane. The energy E associated to a given 
hit is: 



E 



CW 



Q(t) 



(2) 



where C is the calibration factor; W is the average energy 
needed for the creation of an electron-ion pair (i.e. 23.6 
eV); R is the electron-ion recombination factor and Q(t) 
is the corrected charge by the electron lifetime from (fj) 
in Sec. mm 

The calibration factor C converts the detector measur- 
ing units (ADC counts) into charge units (fC). To mea- 
sure the overall calibration factor a = ^W- in (J2J), we have 
collected a sample of around 3000 through-going muons 
coming from CC interactions. The deposited charge 
per unit of length (dQ/dx) is precisely measured (see 
Fig. El and fitted to a convoluted Landau-Gaussian func- 
tion to obtain the most probable energy loss for mips in 
terms of detector charge units: (dQ/dx) mi = 216.7±0.1 
ADC counts. The most probable energy loss is a better 
defined quantity than the average value in the case of a 
Landau distribution, moreover it slightly depends on the 
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FIG. 4: Example of the 3D reconstruction of a low-multiplicity CC event. The raw image from collection and induction 
wire planes (on top): Hits and 2D track projections have been identified. Three dimensional view of the reconstructed event 
(bottom picture) embedded in a 3D recreation of the experimental setup described in Sec. ITT1 



energy in the range of the muons being considered. From 
this measurement, a can be obtained assuming that ac- 
cording to the theoretical prediction for muons in the 
range of energies of the beam, the most probable energy 
loss (dE/dx) mip = 1.736 ± 0.002 MeV/cm: 



reason why we take in J2J R m ip (its value for mips). The 
dependence of R with dE/dx can be modeled by Birk's 
law j2(J which has been successfully applied to LAr de- 
tectors in l27ll: 



l imp 



cw 



(dE/dx) 



dQ 
dx 



dE/dx 



imp 



1 + k B dE/dx 



(4) 



(dQ/dx) 



(3) 

The recombination factor R in J5J depends on the ab- 
sorber medium, on the applied electric field and on the 
density of released charge, i.e. on dE/dx. This is the 



= 8.01±0.01xl0~ 3 MeV/ADC 

dE/dx being the pure Bethe-Bloch ionization loss, a the 

energy to charge conversion factor and k B the Birks co- 
efficient accounting for quenching due to recombination. 

The size of the range of the fully contained protons 
collected in the "golden sample" considered in Sec. [V] 
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FIG. 5: (Top) The raw images of a high multiplicity event (z/ M CC DIS) in the collection (left panel) and induction plane (right 
panel). Below, the 3D reconstruction of a high multiplicity event (y^ CC DIS) with eight primary particles. Two of them stop 
in LAr and are recognized as protons. 



offers the opportunity to precisely test the energy loss 
pattern in the active medium. Fig. shows the observed 
charge per unit length (dQ / dx) versus the residual range 
for the protons 2 . Dots are direct measurements from the 
reconstructed hits of the proton tracks. From this dis- 
tribution we estimate the most probable value of dQ/dx 
for each bin of the residual range. Finally, from this data 



The residual range is the actual range minus the range already 
covered at the time of the deposition of a given dQ/dx. 



we fit k B = 0.035 ± 0.001 cm/MeV in Q for the proton 
hypothesis (continuous line). 



Based on the knowledge of the calibration parameters, 
a precise determination of the deposited energy can be 
obtained through direct charge measurements by means 
of Moreover, pure Bethe-Bloch ionization loss can be 
corrected by the detector response in the region of high 
dQ/dx, where quenching effects are more sizable. This 
correction is important for particle identification in LAr, 
which exploits the different energy loss of particles near 
the stopping point. 
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FIG. 6: dQ/dx of crossing muons. The distribution is fitted to a convoluted Landau-Gaussian distribution. 



C. Proton reconstruction 

Proton identification and momentum measurement 
were performed using only information provided by the 
TPC. The discrimination between protons and charged 
pions is performed exploiting their different energy loss 
behavior as a function of the range. For the case of can- 
didates stopping in the fiducial volume of the chamber, 
as the ones of the "golden sample", p/ir^ separation is 
unique (see Fig. |8j [left]): dQ/dx is measured along the 
proton candidate track and is compared with the different 
particle hypothesis. Monte-Carlo studies reveal that al- 
most 100% of protons are identified as such on the basis of 
their dQ/dx shape on the vicinity of the stopping point; 
in addition, the fraction of pions and kaons misidcntificd 
as protons is negligible. 

Once we identify a contained proton, its kinetic energy 
is calculated from range, which only depends on the spa- 
tial reconstruction of the tracks. In this case, the momen- 
tum uncertainty is dominated by the finite pitch of the 
wires (2.54 mm) 3 and the performance of the reconstruc- 
tion tools. These have been tested using a Monte-Carlo 
simulation of the detector. The resolution on the kinetic 
energy (T p ) varies from 3.3% for protons with T p = 50 
MeV up to 1% for T p larger than 200 MeV. The angular 



3 The equivalent pitch in the vertical direction (drift direction) is 
much smaller due to the high sampling rate of the fast ADC 
(360 fim) 



resolution in the collection or induction view depends on 
the number of wires N hit by the particle along its path; 
it is a ~ 0.36%/l2/(2.54 TV 3 / 2 ), corresponding e.g. to 15 
mrad for N=10. 



D. Muon reconstruction 

The kinematic reconstruction of the outgoing muons 
exploits the tracking capability of NOMAD. An event 
triggering the chamber will have at least one penetrating 
track reaching the Tl and T2 trigger scintillators bracket- 
ing the TRD of NOMAD 0. The corresponding track, 
nearly horizontal at the entrance of the NOMAD drift 
chamber volume, is reconstructed with an average mo- 
mentum precision of a p /p ~ 0.05/V^© 0.008p/L 5 / 2 , L 
being the visible range in the volume itself expressed in 
meters and p the particle momentum in GeV. A 10 GeV 
horizontal muon crossing all the chambers (L ~ 5m) is re- 
constructed with a precision of 2.2%. The reconstructed 
particle is traced back to the TPC accounting for the 
magnetic field and the presence of the forward NOMAD 
calorimeter. The latter introduces an additional - in fact, 
dominant - source of uncertainty due to multiple scatter- 
ing (MS) in iron (190 cm for a horizontal muon). For 
small scattering angles (0^1 rad), the MS uncertainty 
on the transverse momentum is independent of p and 
it turns out to be ~ 140 MeV. The correctness of the 
back-tracing procedure has been cross checked comparing 
the direction angles of the particles belonging to the mip 
sample, as measured by the TPC, with the correspond- 
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FIG. 7: dQ/dx as a function of the residual range for protons. Dots are direct measurements from the reconstructed hits of the 
proton tracks, while the points with errors bars are the estimated most probable of dQ/dx for each bin of the residual range. 
Data is fitted with a Bethe-Bloch corrected by the detector response function (Birk's law) in the hypothesis of pion (dashed 
line) and proton (full line) 




residual range (cm) dE/dX (MeV/CITl) 



FIG. 8: (left) dQ/dx as a function of the residual range for the reconstructed proton in the event shown in Fig. [I] continuous 
(dashed) line represents the expected behavior of protons (pions) in the detector, (right) dE/dx distribution for the muon 
reconstructed in Fig. 2] data is fitted to a convoluted Landau- Gaussian distribution. 
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ing quantity from NOMAD 4 . Due to the large amount 
of material between the TPC and the NOMAD trackers 
(more than 11 interaction lengths) and the small dis- 
tance between the chamber and the calorimeter, the tt//i 
misidentification probability is negligible for the present 
study, even accounting for possible tt decays in flight. 



V. ANALYSIS OF QUASI ELASTIC ^ 
INTERACTION 

A. Data taking and event selection 

The data set recorded with the 50L chamber amounts 
to 1.2 x 10 19 protons on target. As mentioned in SeclTTI 
the trigger efficiency was monitored during data taking 
and its integrated value is 97%. Additional losses of 
statistics are due to the TPC (3%) and NOMAD (15%) 
dead time and to detector faults. These contributions 
add up to give an effective lifetime of 75%. Over the 
whole data taking period around 70000 triggers were col- 
lected in the TPC from which 20000 have at least a recon- 
structed muon possibly transversing the fiducial volume 
(see Sec. IIV Dl . From all these a priory CC candi- 
dates, a visual scanning reveals that around half of them 
show a vertex in the fiducial volume, the rest are v inter- 
actions in the surrounding LAr bath (which also give trig- 
ger) plus a contamination of crossing muons due to veto 
inefficiencies. Therefore, we have a collection of around 
10000 CC events from which we have selected a set of 
86 events called the "golden sample" . This set consists of 
events with an identified proton of kinetic energy larger 
than 40 MeV fully contained in the TPC and one muon 
whose direction extrapolated from NOMAD matches the 
outgoing track in the TPC (see Sec II VI) . The distance of 
the primary interaction vertex to any of the TPC walls 
has to be greater than 1 cm. The muon candidate track 
projected onto the wire plane must be longer than 12 
wire pitches. The event is accepted even in the presence 
of other stopping particles, as far as their visible range 
does not exceed the range of a proton of 40 MeV kinetic 
energy. If other tracks than the identified muon leave 
the TPC or at least one converted photon with energy 
greater than 10 MeV is present in the fiducial volume, 
the event is rejected. The tightness of these selections 
makes the "golden sample" a very clean topology for vi- 
sual scanning (Fig.0J). 

The request of containment is very severe since the 
chamber volume is small. For lower T p values, the pro- 
ton range is comparable with the wire pitch and neither 
the proton momentum nor the interaction vertex can be 



reconstructed with due precision. However, for T p > 40 
MeV the 7r ± /p misidentification probability is negligible 
(see Sec. IIV C|) . The "golden sample" contains pure QE 
interactions [y^ n — > /i~ p) plus an intrinsic background 
dominated by resonant production followed by pion ab- 
sorption in the nucleus ( p — ► A ++ /i~ — > p tt + , 
n — > A + /i~ — *■ /i~ p it ). There is also an instrumen- 
tal background due to final state neutral particles (w°'s, 
n's and 7's) that scape undetected. On the other hand, 
the v M n — > A + /i~ — > fi~ n 7r + contamination is neg- 
ligible in this tightly selected sample due to the superb 
tt /p identification capabilities of the LAr TPC. 

B. The Monte-Carlo data sample 

The selection efficiency for QE interactions and their 
intrinsic background have been evaluated using a Monte- 
Carlo sample generated using the FLUKA package [2^, 
l2flj . which offers a full description of nuclear effects in 
neutrino interactions. About 18000 CC QE events 
have been generated simulating the WANF beam for the 
50L chamber geometrical acceptance and final state par- 
ticles were tracked in LAr using the GEANT4 package 
[3^. The energy deposition of all final particles is dig- 
itized emulating the TPC wire read-out, in order to re- 
produce detector resolutions and offline reconstruction 
efficiencies (see Sec. II VL We saw that 16 % of the 
whole sample of generated QE interactions belong to the 
"golden sample" . On the other hand, we have analyzed 
a sample of 10 5 Vn CC deep-inelastic (DIS) and reso- 
nant (RES) events in order to evaluate inefficiencies of 
the vetoing selections for n — > A + /i~ — > fi~ p tt° 
(i.e. the probability to miss both the decay photons of 
the 7T° — * 77 or the e + e _ 7 system in case of tt° Dalitz 
decay) and for events with several neutral particles (neu- 
trons and energetic photons from nuclear interactions) 
which scape from the detector (see Fig. . The search 
for this irreducible background in the DIS plus RES sam- 
ples reveals that 0.14% are "golden-like" events. 

Taking into account the relative weight between QE 
and non-QE events 5 , we find that the total expected con- 
tamination of the "golden sample" is 20% (around 40% 
of them correspond to events with an escaping tt ). The 
Monte-Carlo prediction for escaping tt° can be checked 
using test samples that consist of golden events with one 
(Ni) or two (N2) converted photons pointing to the inter- 
action vertex. Assuming the gamma identification prob- 
ability to be uncorrelated for the two photons, we have 
Ni = 2N(1 - e 7 )e 7 and N 2 = Ne^, N being the (un- 
known) overall rate of p \x~ tt° final states and e 7 the 



4 Due to the high sampling rate of the TPC, the angular reso- 
lution in the y — z plane, i.e. in the vertical plane along the 
nominal beam direction (z-axis), is dominated by the NOMAD 
uncertainty. 



5 Respect to the total number of CC events, 2.3% are QE, 
91.5% are DIS and the remaining 6.2% corresponds to resonances 
events. The simulation of the latest is not done treating each 
resonance individually, an average of them is taken into account 
instead. 
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photon identification efficiency. After the scanning of the 
test samples, e 7 = (N 1 /2N 2 + 1)~ 1 turned out to be (43± 
9) %. Hence, the probability of missing both gammas is 
(32±10) % to be compared with the MC calculation of 
20.4%. 



C. Event rates 

We use Monte-Carlo to estimate the expected number 
of events in the detector. The simulation of the beam 
predicts a flux of 2.37 x 10~ 7 i^CC/em 2 /p 

.o.t. over the 

TPC exposing area. This flux convoluted with the neu- 
trino cross sections and scaled to the fiducial mass of the 
detector gives an event rate of 2.05 x 10 _15 i^CC/p.o.t. 
Now assuming a total exposure of 1.2 x 10 19 p. o.t. and an 
effective lifetime of 75% (see Sec. IV All , the total number 
of i/ p CC is equal to 18450. Out of this 18450 events, 
only muons above a certain momentum threshold and 
within angular acceptance will trigger the detector. Tak- 
ing the real data muons that triggered the system, we 
saw that most of them are above 8 GeV in momentum 
and below 300 mrad in angle. If we apply those cuts to 
the Monte-Carlo samples, we end up with 18450 x 0.023 
(QE fraction) x 0.95 (efficiency of the acceptance cuts 
for QE), equal to 400 quasi-elastic events. For DIS+RES 
events: 18450 x 0.977 (fraction of DIS+RES) x 0.65 
(efficiency of acceptance cuts) gives 11700. In total we 
expect 12100 events to be compared with the 10000 
CC we have after visual scanning of good triggers (see 
Sec.lVAl 

The expected number of "golden events" is obtained 
taking the 16% (golden fraction) of the total QE, and 
adding the corresponding 20% due to background contri- 
bution (see Sec. IV B|) . which finally gives 80 ± 9(stat) ± 
I3(sys) "golden events" to be compared with the 86 we 
observe. The systematic uncertainty is dominated by the 
fraction of QE events (2.3%) and the beam simulation 
(8%). 



D. Analysis of quasi-elastic interactions 

In spite of the limited statistics (86 events), the "golden 
sample" provides information on the basic mechanisms 
that modify the kinematics of neutrino-nucleus interac- 
tions with respect to the corresponding neutrino-nucleon 
process. Nuclear matter perturbs the initial state of the 
interaction through Fermi motion; it also affects the for- 
mation of the asymptotic states through nuclear evapo- 
ration, hadronic re-scattering or hadronic re-absorption. 
Several kinematic variables are only marginally affected 
by nuclear effects; in this case, the corresponding distri- 
butions can be reproduced once the j/-nucleon interac- 
tion is corrected for Fermi motion and Pauli blocking. 
Clearly, purely leptonic variables belong to this category. 
There are also a few hadronic variables whose distribu- 
tion is strongly influenced by the selection cuts but show 



limited sensitivity to nuclear effects. In particular, the 
proton kinetic energy is bounded by the T p > 40 MeV 
cut and the requirement of full containment in the fidu- 
cial volume. This distribution is shown in Fig. 1101 to- 
gether with the transverse momentum of the proton for 
the "golden sample" . Similarly, in Fig.^Jwe show the ki- 
netic energy and the transverse momentum distributions 
of the muon. Figs.HUlandlTTlcan be used as a consistency 
check. They demonstrate that Monte-Carlo reproduces 
the kinematic selection performed during the scanning 
and analysis of the "golden channel" . They also show 
that our selected sample contains a non-negligible con- 
tamination from non-QE events. As already indicated, 
we estimate this contamination to be 20 %. 

Other variables embedding the reconstructed kinemat- 
ics of the protons are sensitive to genuine nuclear effects. 
In particular, we analyzed two of them: the acollinearity 
and the missing transverse momentum of the event. The 
former is defined as 



.4 



Pxp Pxfi H~ Pyp Py^ 
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(5) 



p xp and Pyp being the transverse momentum components 
of the proton and p xfl and p yp the corresponding quan- 
tities for the muon. For a purely QE scattering on a nu- 
cleon, the muon and the proton ought to be back-to back 
in the transverse plane so that the acollinearity is zero. 
Fig. E| (right) shows the acollinearity distribution for the 
"golden sample" and the expectation from simulations. 
For this particular variable, the inclusion of nuclear ef- 
fects in the Monte-Carlo does not show a striking differ- 
ence with respect to the case where no nuclear effects are 
taken into account. The selection cuts (in particular, the 
fact that the proton should be fully contained) , the reso- 
lution NOMAD has for muon reconstruction and the con- 
tamination from non-QE events are the main reasons for 
the appearance of events with large acollinearity values. 
The distortion introduced in the tail of the acollinearity 
distribution by nuclear effects is much smaller than the 
one due to detector effects. Hence this variable is not 
adequate to show how the event kinematics varies in the 
presence of nuclear matter. 

The influence of nuclear effects on the event kinematics 
is best seen when we consider the missing transverse mo- 
mentum p™ lss (see left plot on Fig. I12fl . We observe that 
a naive approach that takes into account Fermi motion 
and Pauli blocking disregarding nuclear effects (dashed 
histograms in Fig. ^| left) does not reproduce the p™ ISS 
data. Once nuclear effects are added, simulated events 
show a good agreement with data. 

In summary, we have seen that nuclear effects are an 
important source of perturbation for the kinematics of 
quasi-elastic neutrino interactions, however the smallness 
of the accumulated statistics does not allow a systematic 
survey of these effects j3l|. Still, these results empiri- 
cally demonstrate the effectiveness of LAr imaging in the 
characterization of low-multiplicity multi-GeV neutrino 
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Collection wires. (128 wfres: 32 cm.) 
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FIG. 9: The raw image of a low multiplicity event in the collection (left) and induction plane (right). The event is reconstructed 
as {v/j, n — > fi~ A + — > p~ p 7r°) with a mip leaving the chamber, an identified stopping proton and a pair of converted photons 
from the 7r° decay. When these photons escape from the chamber, the event is tagged as a "golden event" . 



Proton Kinetic Energy 



Proton Transverse Momentum 



■ 


50L TPC data 




golden QE 


□ 


golden DIS+RES 


□ 


golden QE+DIS+RES 




0.05 0.1 0.15 0.2 0.25 0.3 

T (GeV) 



(5 

i- 35 



> 
LLI 



■ 50L TPC data 

golden QE 

^ golden DIS+RES 
j | golden QE+DIS+RES 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Pt (GeV) 



FIG. 10: Distribution of the proton kinetic energy for the "golden sample" (squared marks) . Dotted histogram represents the 
expectation from simulated Monte-Carlo QE events, while the hatched one represents the simulated background from non-QE 
events. Both contributions are summed in the filled histogram and normalized to the golden sample data. 



interactions. In particular, a good description of trans- 
verse variables is fundamental since they are key to a 
kinematic-based approach for — > Vrivg) oscillation 
appearance |32|. 



VI. CONCLUSIONS 

We discussed the first exposure of a LAr TPC to a 
multi-GeV neutrino beam. The data collected at the 
CERN WANF exploited the 3D reconstruction capabili- 
ties of the LAr technology with low- multiplicity events. 
The instrumental performances (size and stability of the 
electron lifetime, reproducibility of the charge response 



for mip and highly ionizing particles) and current off-line 
reconstruction tools show that LAr TPC allows for an ex- 
cellent particle identification and a precise measurement 
of the event kinematics. In spite of the limited size of the 
detector, nuclear effects beyond Fermi motion and Pauli 
blocking (e.g. hadron re-scattering or pion absorption) 
have been observed as perturbations of the quasi-elastic 
Vfj, CC interaction kinematics. Variables in the transverse 
plane, which are the most sensitive to nuclear effects and 
key for neutrino oscillation searches, are well reproduced. 
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FIG. 11: Distribution of the muon energy and the muon transverse momentum for the "golden sample". Dotted histogram 
represents the expectation from simulated Monte-Carlo QE events, while the hatched one represents the simulated background 
from non-QE events. Both contributions are summed in the filled histogram and normalized to the golden sample data. 




FIG. 12: Distribution of the missing transverse momentum (left) and acollinearity (right) for the "golden sample". Dotted 
histogram represents the expectation from simulated Monte-Carlo QE events, while the hatched one represents the simulated 
background from non-QE events. Both contributions are summed in the filled histogram and normalized to the golden sample 
data. The dashed histograms show the expected distributions in case no nuclear effects are taken into account. 
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